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Abstract 

New aspects of spectral fluctuations of (quantum) chaotic and diffusive sys- 
tems are considered, namely autocorrelations of the spacing between consecu- 
tive levels or spacing autocovariances. They can be viewed as a discretized two 
point correlation function. Their behavior results from two different contribu- 
tions. One corresponds to (universal) random matrix eigenvalue fluctuations, 
the other to diffusive or chaotic characteristics of the corresponding classi- 
cal motion. A closed formula expressing spacing autocovariances in terms of 
classical dynamical zeta functions, including the Perron-Frobenius operator, 
is derived. It leads to a simple interpretation in terms of classical resonances. 
The theory is applied to zeros of the Riemann zeta function. A striking corre- 
spondence between the associated classical dynamical zeta functions and the 
Riemann zeta itself is found. This induces a resurgence phenomenon where 
the lowest Riemann zeros appear replicated an infinite number of times as res- 
onances and sub-resonances in the spacing autocovariances. The theoretical 
results are confirmed by existing "data". The present work further extends 
the already well known semiclassical interpretation of properties of Riemann 
zeros. 
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I. INTRODUCTION 



Spectral fluctuations of classically chaotic or diffusive quantum systems have been ex- 
tensively studied in the past. One of the main outcomes of these investigations has been to 
establish energy (or time) scales for which universal properties hold and which are repro- 
duced by random matrix theories (RMT), to be distinguished from those which are system 
dependent and in many cases can be treated by semiclassical theories a la Gutzwiller#. 
Apart from the one-point function or level density p which gives the global behavior and 
sets the main scale (its average p defines the mean spacing D = p" 1 , or its conjugate variable 
defines the Heisenberg time T# = hp, with h the Planck constant), much emphasis has been 
put on the two-point function R%{e) or density-density correlation function or on quantities 
therefrom derived (R2(e)de is proportional to the probability of finding two levels separated 
by a distance in the interval [e, e + de]). One important exception to this is constituted by 
the nearest neighbor spacing distribution p(s) {p{s) is not a two-point function). It has been 
thoroughly investigated and its exact short and long range behavior have simple analytical 
forms for the Gaussian Ensembles. The short range behavior of i?2(e) at the scale of one 
mean level spacing is the same as the one of p(s), whereas the information contained in the 
medium and long range part of i?2(e) (medium and long range correlations) is not contained 
in p(s). It is rather well established that semiclassical theories are particularly adapted 
to describe medium and long range correlations whereas their ability in dealing with short 
range properties is dubious. 

Studies of autocovariances C (n) of the spacing between consecutive levels, namely quan- 
tities related to the average value of s m s m+n , where s m = t m — £ TO _i, are scarce (t m denotes 
an eigenvalue, with . . . < t m < t m+ i <...). As will become clear in section II. A, spacing 
autocovariances C(n) can be viewed as a discrete version of the (continuous) two-point cor- 
relation function (cf Eq.( |2.6| )). This is related to the fact that the basic building block is a 
discrete object, a spacing between consecutive levels. The discretization on a scale of the 
mean level spacing induces a "smoothing" procedure, and the structures existing up to that 
scale are strongly suppressed. As a consequence the spacing autocovariances are particularly 
suited to be described by semiclassical approaches, which are usually very difficult to con- 
trol on the mean spacing scale. In terms of the conjugate variable (time), the discretization 
amounts to introduce a cut-off at the Heisenberg time T#. 

In contrast to R2{c), exact analytical random matrix expressions for the spacing auto- 
covariances C(n), whose study is one of the main purposes of the present paper, are not 
known (this explains why they are usually not extracted from data). In order to proceed 
use will be made of a simple relation connecting spacing variances of distant eigenvalues to 
the number variance (Eq.(|2.5|)). As discussed in Ref.(Q), this very general relation applies 
to eigenvalues of Gaussian Ensembles and also to spectra of chaotic systems, even beyond 
the universal regime, as well as to spectra of integrable systems (but does not hold for Pois- 
sonian spectra). It is therefore associated to the concept of spectral rigidity, characteristic 
of Gaussian Ensembles of random matrices and of the non-universal regime of quantized 
Hamiltonian systems. 

This paper, which may be considered as a continuation of Ref.(|2|), is organized as fol- 
lows. The next section starts with the general setting, definitions and relations which will 
be subsequently used (section II. A). In particular, we establish an approximate (but accu- 
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rate) relation between spacing autocovariances and the counting function variance (number 
variance). When applied to eigenvalues of random matrices the resulting expressions are 
consistent with an ansatz introduced in Ref.(|2]), though the precise value of constants differ 
because of slight differences in both treatments. Quantum systems governed by classically 
diffusive motion (section II. B) and by chaotic dynamics (section II. C) are subsequently stud- 
ied. In the latter case the main tool used is a semiclassical approach. Because the quantities 
computed are less sensitive to small energy scales, we find that the simplest form of the 
theory, i.e. the "diagonal approximation", provides a very accurate global description of 
C(n). The main result of section II. C is an expression of the autocovariances in terms of 
two different classical dynamical zeta functions, one of them directly related to the spectrum 
of the Perron- Frobenius operator. The structure of C (n) can thus be interpreted as a super- 
position of contributions of classical singular points (or resonances). Section III has initially 
been conceived as an illustration of what precedes. One studies properties related to the 
zeros of the Riemann ^-function by interpreting them as eigenenergies of an hypothetical 
quantized classically chaotic system, as suggested by Berry exploiting a beautiful analogy 
with the Gutzwiller trace formula®. Expressions for the spacing autocovariances are worked 
out and compared to "data" provided by the extensive computations of Odlyzko®. Fea- 
tures encoded in spacing autocorrelations are exhibited. Most of them are new and some 
unsuspected. That is why we think that the content of the section has interest in its own 
and goes beyond its initial illustrative purpose. 

II. SPACING AUTOCOVARIANCE 
A. General Considerations 

Consider an ordered sequence of energy levels ti, % = 1, 2, . . ., having an average density 
of states p{t) . We call them "energies" , but the set can be an arbitrary sequence of points on 
the real line not necessarily related to a quantum mechanical spectrum. From this sequence 
a new stationary one, having mean level spacing equal to one, is constructed by the unfolding 
or rectifying procedure®'® Xi = /** p(t)dt. The sequence x» is located around energy t in a 
window of size At, with At << t. 

Our purpose is to study properties related to distances between two consecutive eigen- 
values, denoted s m = x m — x m -\- More precisely, we study the autocovariances C{n) of two 
spacings between consecutive eigenvalues located n levels apart 

C(n) = \{s m — (s m ))(s m+n — (s m+n )) = (s m s m+n ) — 1 = I[n) — 1 , (2-1) 

where I(n) are the spacing autocorrelations. The brackets denote an average over different 
locations m of the reference spacing or, more generally, an energy average over the spectrum. 
Eq.( [2.1|) is the simplest quantity containing information on how the different spacings are 
correlated that can be considered (correlations among spacings of non-consecutive levels 
can be expressed in terms of C(n) and do not contain new information). For uncorrelated 
spacings C(n) = (s m )(s m+n ) -1 = 0. 

Consider the length 5* of an interval made of n consecutive spacings, 5* = J27=i s i ■ The 
variance of 5* can be written as 
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n-1 



a 2 (n) = {S 2 )-{S) 2 = n ((s 2 )-n)+2 ^(n-j)I(j), n>2 



(2.2) 



and cr 2 (l) = (s 2 ) — 1. Inversely, Eq. 
the variance of S, 



allows to express the autocovariances in terms of 



C{n) 



1 r 

2 
1 



ff 2 (n + l)-2ff 2 (n)+ a 2 (n-l) 



n > 2 



C(l) = -[a 2 (2)-2a 2 (i; 



(2.3) 
(2.4) 



For eigenvalues of Gaussian Ensembles of random matrices neither a 2 (n) nor C(n) are known 
in closed form. Our aim is to show that an approximate (but accurate) expression for the 
autocovariances of spacings may be found exploiting some properties of a 2 



,n 



For that purpose we consider the number variance S 2 (L) closely related to a 2 (n ) 



and Eq.(|J). 
.2^ Itig 

defined as the variance of the random variable which counts the number of levels contained 
in an interval of length L located randomly in the sequence. For eigenvalues of Gaussian 
Ensembles of random matrices, S 2 (L) evaluated at integer values of its argument is related 
to the variance of the interval length throug h§ 



E 2 (L 



n) 



a 2 (n) 



1/6 



(2.5) 



This relation is valid, in principle, for large n. However, numerical calculations as well as 
analytical estimates© indicate that Eq. fl2.5| ) holds also for small values of n with an error 
~ 0.01/n 2 . Eq. ([2.5|) was used in Ref.(0) to determine the coefficients of an ansatz for C{n) in 
terms of a sum of inverse even powers of n. In this section we follow an alternative path. We 
make no ansatz, but use Eq.( |2.5| ) to derive similar random matrix theory (RMT) expressions 
for the autocovariances C(n). The next sections will be devoted to the incorporation of the 
semiclassical and diffusive contributions to the autocovariances. 
Assuming that Eq. ( |2.5|) holds, Eqs. ( |2.3| )- (|2T4T ) can be written 



C(n) 
C(l) 



E 2 (n + 1) - 2 S 2 (n) + S 2 (n - 1) 



n > 2 



E 2 (2) -2 E 2 (l) + 



(2.6) 
(2.7) 



Due to the structure of Eq.( |2.3| ) the exact value of the constant in the r.h.s. of Eq. ( |2.5| ) 
is immaterial for n > 2. Notice that Eqs. (|2.5| ) and (|2.6j ) connect quantities related to the 
exact location of energy levels, like a 2 (n) or C(n), to quantities which are not, like S 2 , and 
for which explicit results exist for chaotic and diffusive systems®'®. 

Using the approximation ( |2.6| ) the autocovariances C(n) are expressed in terms of the 
(discrete) curvature of S 2 evaluated at integers n. An analogous expression (up to a sign) 
relates for L > the two-level cluster function® Y 2 (L) to the (continuous) curvature of 
S 2 (L), 



l d 2 E 2 

2 ~mJ 



(2.8) 



4 



One therefore expects the properties of C(n) to be closely related to the ones of Y 2 . One 
aim of the present paper is to identify some of their differences. 

Let us recall that the Fourier transform of Y 2 defines the form factor K(r), 

POO 

Y 2 (L) = 2 / dr(l - KM) cos(27rLr) , (2.9) 
Jo 

in terms of which the number variance reads 

J?(L) =* r ir K(t)**¥1 . (2.10) 

7T Z JO T z 

t is a rescaled time in units of the Heisenberg time, r = T /hp and, correspondingly, energies 
are being measured in units of mean spacing p _1 . 

As is well known, the form factor K(t) plays an important role in the theory of quantum 
dynamical systems. For disordered and chaotic systems, it contains basic information related 
to short time specificities. In the opposite limit, namely the ergodic (long time) regime, its 
behavior is universal and coincides with RMT. Accordingly, in order to exhibit the different 
relevant time scales, it is convenient to write 



+ £nu- (2-11) 



i v*2 

J rm 



E^ m is the random matrix result for the number variance (cf Eq. fl2.15p below). The non- 
universal contribution will be considered later on for disordered as well as for chaotic 
systems. From Eq. fl2.6|) , and in analogy with Eq.( |2.11| ) we write 

C = C rm + C nu . (2.12) 



Using the exact form of it is possible to obtain, through Eq. (^6|) , a general expression 
for C rm . However, we rather prefer to use an asymptotic expansion of the number variance 
for large values of n. In particular, this makes easier the comparison of the present approach 
with the results of Ref.(^). For large values of L, the leading order behavior of is® 

£ 2 (L) = -^logL + 0(l), (2.13) 

where j3 =1,2 and 4 denote the three types of Gaussian Ensembles, orthogonal, unitary 
and symplectic, respectively. Eq.( [2.13| ), together with Eq.(|2.6|) leads to 

C,W = ^.og(l-l) , (2.14) 

valid for n > 2. For a more accurate evaluation of C(n) we include more terms in the 
asymptotic expansion of the number variance. For example, for (3 = 2 we use 



£ 2 (£) = - 2 



cos(2vtL) 4sin(27rL) „„. , 4n 
log(2,L) + 7 + 1 - -pLi - + OWL*) 



(2.15) 



An analogous expansion can be used for /3 = 1. Substituting these expansions in fl2T6|) one 
gets, for (3 = 1 and 2 
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Crm{n) = CUn) + 0{l/n*) , (2.16) 



with \p and ctp given by 



Ai = 


3 

2tt 2 ' 


= 


15 
6tt 2 


105 

6tt 4 


A 2 = 


3 


a 2 = 


15 


135 


2tt 2 ' 


67T 2 


6tt 4 



(2.17) 
(2.18) 

For (3 = 4 a theorem relating the statistical properties of the Gaussian Ensemble with f3 = 1 
to those of f3 — 4 can be exploited®. It implies that the autocovariances are related as 
follows (the notation is obvious), 

C rm ,4(n) = ~ C rmA (2n) + i C rin ,i(2n - 1) + J C ro ,i(27i + 1) . (2.19) 



Eq. (|2.19|) enables to determine the coefficients in fl2.16|) for (3 



A4 = Ai/4 , «4 = (ai + 5 Ai)/16 . 

Aside from the logarithmic leading order term, Eq. fl2.15|) contains sub-dominant oscil- 
latory corrections. The latter modify the simple estimate Eq.( |2.14|) by adding smooth cor- 



rection terms of order 0(l/n 4 ) and higher according to Eq. (|2.16 ). Contrary to this, when 



computing the cluster function from Eq. Q2.8|) these sub-dominant terms of £ 2 (£) give rise 
to an oscillatory behavior of Y2 on the scale of the mean level spacing. For example, for 
(3 = 2 we have 

y(t\- 1 cos ( 27rL ) room 
Y2{L) ~ 2^ ' ~2^IT ■ (2 ' 20) 

In contrast, there are no oscillations in C(n) because S 2 is evaluated at integer values of its 
argument, and the oscillatory functions in Eq.( [2.15 ) are either zero or one. Similar results 



are also found for the other symmetry classes. This makes an important difference between 
the RMT behavior of Y2 and that of C(n), with a suppression of the structures of Y-x on the 
scale of the mean level spacing or below. 

Expanding the logarithm in C® m (n) one obtains from Eq. fl2.16Q a series in inverse even 
powers of n, which is consistent with the ansatz of Ref.(|2p. However the coefficients in this 
expansion slightly differ from those obtained in Ref.(fJ). This difference is due to the fact that 
terms of 0(l/n 2 ) have been ignored in Eq.( |2.5| ). Despite its asymptotic character, Eq.( |2.16 ) 



reproduces with good precision the autocovariances of spacings of random matrices down to 
values n ~ 2. 

Let us finally mention that the sum-rule 

°° <T 2 (1) 

£CX?) + ^P = (2.21) 

3=1 

valid for Gaussian Ensembles©, is violated in the present approach due to inaccuracies of 



the low (j = 1 and 2) terms (C(l) has to be computed from Eq.(277)). This is in contrast 
with what is done in Ref.(^j) where, by construction, Eq.( |2.2lD is fulfilled. Notice however 
that the procedure adopted here is much better adapted for a semiclassical description and 
that this apparent deficiency is of no consequence for what follows. 
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B. Diffusive systems 



We now consider the motion of a particle confined to a spatial region of typical size R 
whose classical dynamics is diffusive. The classical motion is characterized by the diffusion 
constant D. The elastic mean free time between elastic collisions is denoted r e = T £ /Th 
(we again measure the time in units of Heisenberg time Tjj = hp). Quantum mechanically 
the relevant diffusive time scale is the Thouless time, i.e. the typical time it takes to the 
system to cross the sample, r c = R 2 /(Dhp) (l/r c = g is the dimensionless conductance). 
For a diffusive system r e « r c « 1. 

Ignoring the specificities of the dynamics for extremely short times r <r e @ , for diffusive 
systems the form factor has, as in the chaotic case, two distinct regimes. For times r << r c 
the dynamics is diffusive (non-universal) and the form factor is given by® 

Knu(r) = l(^-) d/2 r 1 ~ d/2 , (2.22) 



with d the space dimension. On the other hand, for r >> r c the ergodic regime is reached 
and we recover a universal behavior for K(r) described by RMT^. We do not intend to 
describe in detail the transition between these two regimes, and simply assume K{r) = 
K TTa (r) + K nu (r) (consistently with Eq.( |2.11D ). 

Inserting K rm (r) in ( |2.10| ) we again get £ 2 m as in the previous subsection. On the other 
hand, from the diffusive form factor ( [2. 22 ) we obtain £ 2 



nu' 



Y 2 m 4 fLr c \ d/2 [°°, sm 2 (nr) 

The integral gives a constant ad = tt,tt 2 /2 and 47r 2 /3 for d = 1,2 and 3, respectively. 
Eq.( [2.23P provides the well known® asymptotic expression of the number variance for 
diffusive systems, £ 2 oc (Lr c ) d ^ 2 . 

Computing the discrete curvature (|2.6p we have 



C — C rm + Cdif . (2.24) 
C rm is the random matrix autocovariance given by (|2.16|) . C^a is computed from ( 2.23 ), 



Notice that for d = 2 the contribution vanishes, and that it is positive for d = 3 and negative 
for d — 1. Though this is a decreasing function of n, its relative importance with respect to 
C rm increases since for large separations it vanishes as n~^ 2 ~ d ^ 2 \ which is slower than the n~ 2 
decay of C rm . It therefore dominates the tail of the autocovariance function in any dimension 
(except d = 2). It represents a finite conductance correction to the RMT behavior, since as 
mentioned before g = t~ x . For d = 3, C(n) is negative for small n (spacings separated by 
small distances are anti-correlated) but tends to zero from above for large distances (positive 
correlations between large distant spacings (see Fig. 1)). For d — 1 C(n) is negative and 
spacings are anti-correlated for any n. 
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C. Chaotic systems 



For ballistic fully chaotic systems K(t) has two distinct regimes®. For times r « r min , 
with r min = T min /hp the rescaled period of the shortest periodic orbit, the behavior of 
K(t) is non-universal, with peaks located at the periods of the shortest periodic orbits. 
For longer times the average behavior of K(t) is universal and consistent with RMT. The 
transition between these two regimes is described semiclassically by the Hannay-Ozorio de 
Almeida sum rule©, and occurs around a cut-off time satisfying r min « r* << 1. 
Correspondingly, for energies n^l/r*, S 2 (n) agrees with the random matrix prediction, 
while around n ~ 1/r* there is a transition towards a non-universal regime, where oscillations 
associated to short periodic orbits around a saturation value occur. These are described® 
by the non-universal term in ( p. 11 ) 



with 



y 2 _i_ y 2 
^po ^ ^* > 



/?7T 2 



E 

rT„<r* 



det(Ml 



sm TcnrTr, 



(2.26) 



and 



[Ci(27rnT*) — log(27rnr*) — 7] 



^2.27) 



In Eq. 

1,2,3, 



the sum is over all the classical periodic orbits p of the system, with r = 
the number of repetitions, t p = T p /Th their normalized period and M p the mo- 
nodromy (or stability) matrix. The sum runs over orbits whose period (including repetitions) 
rr p is smaller than r*. 

From (g^BJ) and fl227|) we get (cf fl2TTj )) 



C C n 



Crm Cr- 



with 



/3tt 2 



E 

TTp <T* 



sm 7rrr, 



det(M; - 1] 



cos(27rnrr p 



(2.28) 



and 



C m (n) 



(3^ 



Ci(27r(n + l)r*) - 2Ci(27mr*) + Ci(27r(n - l)r*) - log 1 - 



^2.29) 



The sum C po + C* accounts for the difference between the short time dynamics of a 
real physical system and a pure RMT description. Indeed is, with opposite sign, the 
RMT contribution for times < r < t* where a linear (diagonal) approximation of the 
form factor has been used. On the other hand, due to the Hannay-Ozorio de Almeida sum 
rule, for times r > the contribution of C po coincides with — C*. This means that we may 



extend r* to infinity in C po and C*, without affecting their sum. But from Eq.( |2.29D we have 
linx^ooC* = -(/37r 2 )- 1 log(l - 1/n 2 ). It follows that 
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sin 2 (7rrr T) ) 1 , / I 



n- 



In contrast to Eq. (|2.28|) , there is no restriction imposed in the latter equation and the sum 
extends over all periodic orbits p and their repetitions. Notice that the second term in the 
r.h.s. of Eg. ( |2.30| ) cancels the term C° m of the universal contribution C rm in Eq.( |2.16|) . The 
autocovariance therefore takes the form (for n > 2) 

. . 4 ^ sin 2 (7rrr„) . . 1 ( \a oin\ m , , S s 

C(n) = W £ >-|det(W-l)l + W U + + 0(V " 1 ■ <2 - 31) 



Eq. fl2.31|) indicates that the essential features of the autocovariances are described by the 



"diagonal" approximation given by the sum over the periodic orbits, which includes the 
universal term C® m given by Eq. fl2.14|) as well as the non-universal fluctuations due to 



shorter orbitsQ What is left out of the sum are the smooth remaining RMT terms of order 
1/n 4 and higher appearing in Eq. fl2.31|) , which are small compared to the leading order 
behavior of C° m ~ — (/?7r 2 n 2 ) — 1 because of the size of the constants and otp (they are of 
order 10 -1 ), combined to the fact that n > 2 in Eq. ( p.31|) . 

In summary, to a very good approximation the autocovariances C(n) are given by the 
first term in Eq. ( |2.31| ) 

an)^E r2| r t 2 ff 1) | OOs(2^ p ). ,,32) 

In this sum, the long, exponentially numerous periodic orbits produce a coherent contribu- 
tion given by Eq. ( j2.14|) , which is significant for low values of n and is independent of the 
particular position or height t of the window considered to compute C{n) (i.e., is independent 
of Th). On the other hand, the amplitude of the contribution of the short orbits scales like 
Tp oc Tft 2 . Therefore, in the extreme semiclassical limit Th — > oo, only the universal term 
contributes. However, at a finite Th the contribution of short trajectories must be taken into 
account, and is particularly important for large values of n where the universal contribution 
is negligible. The distribution of the t p , their arithmetic and commensurability properties 
for low values of p, as well as the stability of the orbits control the form of C(n) for large val- 
ues of n. In the simplest (and probably general) case where only a few non-correlated short 
orbits significantly contribute to C(n), the sum ( p.32[ ) leads to strong long range correlations 
between the spacings s m . The exact form depends on the interference pattern between the 
orbits, with a rough estimate of the amplitude given by 4r^ in //3| det(M min — 1)|. However, 
coherent interference contributions between several orbits and its repetitions are of course 
possible, and larger fluctuations cannot be excluded. 



2 From a numerical point of view, instead of computing the sum in Eq.( 2.3l| ) over all the periodic 
orbits and its repetitions it is more convenient to truncate the sum at some rr p < r#, where t* is 
sufficiently large to ensure a RMT behavior, and add — C* (with C* given by ( p.29| )) to account 
for the remaining part of the sum. 
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In order to clarify this point and to acquire further insight on the structure of C(n) we 
rewrite the sum ( 2.32Q in a different manner. Our purpose is to express the autocovariance 
in terms of classical zeta functions, a procedure that leads to a more transparent analysis. 
For that end, we first expand the sine function in a power series to get 

1 oo 2(fc-l)o2fc+l 2(fc-l) 2fc 

C(n) = T;Y(-l) k+1 , ,7, Y tt— — rr cos(2vrnrr p ) . 

1 ' Pti ( 2k V- £|det(Af;-l)| 1 pJ 

It is convenient to split these sums in two parts and corresponding, respectively, 
to the two terms on the r.h.s. of the decomposition of the 2fc-th power of the period in the 
numerator 

rf = r rf - (r - 1) rf . (2.33) 
The contribution of the first term to C(n) may be written 



(3tt 2 (2k)\ dn 2k j£ r | det (M£ - 1) | 
or, alternatively 

( n ) = —Y -—Re^— log TT exp , — . (2.34) 

fin 2 j^i (2fc)! <9n 2fc P V r | det(M; - 1)|/ v 7 

The argument of the log coincides with a well known classical dynamical zeta function©'© 

/ srT \ 

z(s) = n n c -p - ,JLr tti . ( 2 - 35 ) 



p r=l 

and takes the form 



r |det(M 



C"»W = ^i; iMJ Re^lo g Z(WM. (2.36) 

This expression is similar, for k = 1, to the relation found® between the diagonal part of 
the spectral two-point correlation function and Z(s) (see also the concluding remarks). 
The contribution to C(n) coming from the second term in the r.h.s. of Eq. (|2.33|) , denoted 
can be obtained following similar steps to those leading to Eq.( [2.36D . We obtain 

O OO 1 Q2fc 

^('^-^E^Re^logFOVM, (2-37) 
with F another dynamical zeta function 

™=nn«*(-?feg^) . (2,8) 

By this procedure, the two terms C^\n) and C^(n) are expressed without any approxi- 
mation in terms of classical dynamical zeta functions. Eq.( [2.32j ) can thus be written in the 
compact form 
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f3Tr 2 ]^(2k)r"dn 2k 



Q2k 

Re - „, log 



Z (in /hp) 
F (in /hp) 



(2.39) 



By using the fact that for chaotic systems the periodic orbits are unstable and that r > 2 
in the expression of F(s), it is reasonable to assume | det(Mp — 1)| 3> 1. This leads to the 
approximate expression 



F( S )«nn 1 



e srT p 



p r=2 



det (Mr - 1) 



(2.40) 



Equation (|2.39| ), together with Eqs. (|2.35|) , (|2.38|) or (|2.40|) , is one of the main results of 
this paper. The net outcome has been a transcription of the sum over periodic orbits by 
an expression in terms of classical dynamical zeta functions. The first one, Z(s), is well 
known®'©. Its complex zeros determine the spectrum 7^ of the Perron- Frobenius evolu- 
tion operator of a classically chaotic system. Physically, the characteristic times (Re7 /J ) _1 
determine the relaxation of the time evolution of classical initial clouds. The second zeta 
function, F(s), is formally very similar to Z(s) but has, to our knowledge, not been studied 
so far. 

Assuming Z(s) and F(s) are entire functions that have no singularities other than poles, 
the beauty of Eq.( |2.39|) is fully revealed by considering the analytic structure of the classical 
zeta functions, i.e. their singular points (zeros and poles) also called resonances^. To simplify 
the discussion, consider first the zeros 7^ of Z(s) (other singular points are briefly discussed 
below). Factorizing Z(s) in terms of the 7^, Z(s) oc Y\^(s — 7^), and substituting this in 
Eq.( |2.39| ), it follows that the contribution to C(n) of the zeros of Z(s) is 



C(n) 



1 



/3tt 



[n + in. 



(2.41) 



where the n M are complex numbers determined by the (inverse) of the classical times mea- 
sured in units of Heisenberg time 

?V = hp 7m • 

The structure of C(n) as a function of the real parameter n may therefore be viewed as 
resulting from the contributions of classical resonances located at points n M in the complex 
plane. These are system dependent, except 70 = n = which always exists for ergodic 
systems and brings in the universal contribution Eq. (|2.14|) . The remaining resonances 7 M , 
p > 1 contribute to the non- universal oscillations of C(n). Because they are located at 
a finite position in the complex plane (different from the origin), in the semiclassical limit 
Tu — > 00 all the non-universal resonances are pushed towards infinity (|n M | — > 00 for p 7^ 0), 
thus leaving the universal term no = as the only remaining contribution. When Tjj is finite, 
each of the terms in ( p. 41 ) produces a peak in C(n) whose shape is given by 



3 By an abuse of language, from now on we refer to poles and zeros as singular points. This 
terminology is justified because, up to a sign, from an analytic function point of view they both 
play an equivalent role in Eq.( 2.39| ). 



11 



fix 



2pn- 



log 



1 + 



2n 



fiR 



in 



n^if + 1" 



It R 



+ (n - 7fyj) 2 p 



(2.42) 



where n^R and n^j are the real and imaginary part of n M , respectively. The peak is centered 
at n « n^i, with a height H = (2/?7r 2 ) -1 log[l + 2(n 2 R + l)/n* R ] and width W = 2?^. Away 
from the maximum, f ^ has a small negative tail that tends to zero as — [(37c 2 (n — n^/) 2 ] -1 . 
The position of the peak is therefore controlled by the imaginary part n^i, whereas the real 
part n^n determines its height and width. 

The other singular points of Z(s) and F(s) contribute with similar terms as those of 
Eqs.( |2.41| ) and (|2.42j) . The only difference is a possible overall sign with respect to the zeros 
of Z(s). Due to the structure of Eq. (|2.39|) , the poles of F(s) contribute with the same sign 
as the zeros of Z(s), whereas the zeros of F(s) and the poles of Z(s) have the opposite sign. 
The final expression of C (n) is therefore given by a superposition of the contributions of all 
the singular points n sp of both functions, taking into account their appropriate sign, 



C(n) 



(3n 2 



^sign(sp) log 



sp 



1 - 



[n + in sp ) 2 



(2.43) 



The larger contributions will come from resonances with the smallest real part. 



III. APPLICATION TO THE RIEMANN ZEROS 

We illustrate here the results of section II. C by considering the complex zeros of the 
Riemann zeta function ((s) = J2n n ~ s - This function is becoming an important model in 
the theory of quantized classically chaotic Hamiltonian systems. The connexion with a 
dynamical system arises when the complex zeros of ({s), located by the Riemann hypothesis 
on the line s = 1/2 + it, are interpreted as the eigenvalues of an hypothetical chaotic 
system®-©. 

At present there are many evidences supporting the idea that asymptotically (i.e., for 
large t) the local statistical properties of the Riemann zeros are described by the GUE 
distribution of random matrices®'®, although more precise statements concerning the 
appropriate matrix ensemble have recently been made©. In his extensive numerical studies 
of the statistical properties of the zeros, Odlyzko® observed strong long range correlations 
between the spacings s m . He also found in the Fourier transform of the autocovariances 
clear evidence of the contribution of small prime numbers. Here we make a quantitative 
comparison between his computations and the analytical results of the previous section. 

For that purpose we adapt the equations of section II. C to the Riemann zeros by ex- 
ploiting the analogy between a classical dynamical system and the Riemann zeta function. 
When comparing the Gutzwiller trace formula for the spectral density of a chaotic system 
on the one hand and the density of zeros on the critical line on the other (here interpreted 
as quantum eigenvalues), it appears that the periodic orbits of the hypothetical "Riemann 
dynamics" are labeled by the prime numbers p — 2, 3, 5, . . .. Many properties of these peri- 
odic orbits are known, and Table I is a list of the more relevant ones. In the same table we 
have also specified other quantities and correspondences, like the (leading order) behavior 
of the density of zeros at a height t along the critical axis. 
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With these definitions in hand, the autocovariance (|2.32|) reads 



C(n) = ^ £ Sin2( 2 7r 7 p) cos(27r»rr p ) . (3.1) 

7T p^ r T p 

The sum now runs over the prime numbers p, and r p = \ogp/ log(t/27r) as indicated in 
Table I. The parameter t sets the height on the critical line where the correlations are 
computed. In Figs. 2 and 3 the theoretical prediction (|3.1| ) is compared with numerical values 
("data") computed by Odlyzko. The agreement is good, although some small deviations 
remain. Important differences with respect to the RMT result are observed already at 
n ~ 3. For small values of n (Fig. 2) the typical amplitude of the non-universal fluctuations 
is of order 10~ 3 . For larger values of n (Fig. 3), the amplitude of the oscillations is of order 
10 -2 , confirming the existence of strong long range correlations. Indeed, assuming that 
the s m are uncorrelated variables, the statistical fluctuations due to finite size effects are 
estimated® to be of order 10 -4 . Besides their amplitude, Fig. 3 shows the existence of sign 
correlations between several consecutive points. 

It is instructive to look at the structure of C(n) on larger scales (Fig. 4). C(n) is computed 
from Eq.(|3.1|), which follows very closely the numerical results of Odlyzko. For relatively 



small separations 50 < n < 500, C(n) presents a remarkable behavior, with a dominant 
small positive correlation "punctuated" by large peaks of anti-correlation (made of 3 to 
4 spacings). Qualitatively, the large peaks are produced by the constructive interference 
of a relatively small number of short periodic orbits. Indeed, we find that the qualitative 
features of the main resonances in Fig. 4(a) are reproduced by summing in Eq. ( |3.1| ) over 
approximately 12 prime numbers only (without repetitions). As n increases the isolated 
peaks tend to disappear, as well as the asymmetry of the plot with respect to the horizontal 
axis (cf part (b) and (c) of Fig. 4). As we shall now see, this peculiar structure finds a 
very simple explanation within the general framework based on the classical zeta functions 
presented in section II. C. 

The first step required to write C(n) as a resonance formula corresponding to Eqs.( |2.35D 



and (|2.43|) is to identify the classical dynamical zeta functions associated to the "Riemann 
dynamics". This is achieved by computing Z(s) and F(s) from Eqs.( |2.35D and ( |2.40| ) using 



the periodic orbit correspondences of Table I (we use here for F(s) the approximate equation 
instead of the exact one). For the former one, this procedure leads to 



no 



p- 1 =C'(l-»), (3-2) 



where £(s) = Yl p (1 — p~ s ) 1 is the Euler product (over the prime numbers p) representation 
of C( s )- Thus, when the critical zeros of the Riemann zeta are interpreted as the eigenvalues 
corresponding to the quantization of some classically chaotic motion, the associated classical 
dynamical zeta function Z(s) turns out to be the inverse of the (translated) Riemann zeta 
itself. Concerning the other dynamical function needed, F(s), using again the periodic orbit 
correspondences we have 



e srT p 



n(i-i^n)iJ=II(i-^) = c- 



rs 
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and therefore, using the approximation ( |2.40| ) we obtain 

F(s) = IJ[C(r - rs)]~ {r - 1)/r2 



(3.3) 



r=2 



This shows the remarkable property that in the Riemann case the (approximate) dynamical 
function F(s) is also expressed in terms of C( s )- The dynamical zeta functions are basic 
tools to study the classical evolution. As we have already emphasized, the Perron-Frobenius 
operator controls the time evolution of classical statistical averages^. We have here com- 
puted two of these classical functions for the "Riemann dynamics" which, for completeness, 
have been added to Table I. 

The resulting expression for the autocovariances written in terms of dynamical zeta 
functions follows from Eqs.( p.39|) , (|3.2|) and (|3.3|) 



C(n) 



-E 



1 



7T 



Q2k 

Re - „, log 



C(l - in/p) 



\ (2k)\ dn 2k 



ni((r -inr/p)} 



-\](r-l)/r 2 



(3.4) 



r=2 



To proceed further we shall exploit the analytic structure of Z(s) and F(s). We simply 
need the following factorization formula© 



C(s) 



e (log 2tt-1-7/2)s 



2( s -i)r(i + s /2)^ 



n 



K 



(3.5) 



which shows explicitly the presence of a pole at s — 1, nontrivial zeros k lying in the critical 
strip < Re k < 1, as well as trivial zeros at the poles s = —2m, m = 1, 2, ... of T(l + s/2). 
Assuming the Riemann hypothesis, namely k = 1/2 + it fl , with t M real and coming in pairs 
symmetrical about the real axis, using Euler's factorization formula for T(l + s/2), and 
deriving term by term in the sums, Eq.( |3.4| ) leads, together with Eq.( |3.5| ), to a resonance- 
type formula 



(S r>1 + 1 - r) 



log 



1 



[n + in 0r ) 



-E lo s 



(n + in^ r ) 2 



(n + in mr ) 2 



. (3.6) 



The outer sum in this equation is over the repetitions. Notice a global change of sign 
(reflecting the structure of Eq. ( |3.4j )) of the r = 1 term with respect to the r > 2 ones. The 
three terms inside the square brackets (of which the last two are sums) originate from the 
pole, the critical zeros and the trivial zeros of ((s), respectively. Due to the structure of 
Eq.(|3.4|) the singular points of ((s) appear in C{n) as resonances located at new, shifted 
positions in the complex plane given by 



n 0r = p (1 - 1/r) 

rv = p[l-(l/2-z^)/r] 



(3.7) 



n,, 



p [1 + 2(m+ l)/r] . 



14 



Notice also, as mentioned in section II. C, that the contribution of the pole has the opposite 
sign with respect to the terms associated to the zeros. 

Consider first the contributions of the singular points when r = 1. Then the resonances 
are located at n 01 = 0, n M i = p (1/2 + i t^), and n m \ = p (2m + 3). The pole n 01 in Eq.( p.6[) 
gives the universal term Eq.( |2.14D . The nontrivial zeros n^i are located on a rescaled and 
shifted "critical" line, and have all the same real part, Ren^i = p/2. This explains the 
very interesting and peculiar structure of the autocovariances: for r = 1 each nontrivial 
zero of £(s) contributes to C(n) with a negative correlation peak whose shape is given by 
Eq.( |2.42| ), centered at n — Imn M i = p t M (= 55.1, 81.9, 97.4, . . . for p — 1, 2, 3, . . ., with p = 
3.89533 corresponding to the window used by Odlyzko in his numerical computations). These 
negative peaks are clearly visible in Figs. |] and [|. The height H and width W of the peaks is 
almost constant as observed in part (a) of Fig. |] because the complex resonances are all at the 
same distance p/2. From section II.C we have H = (47T 2 )- 1 log[l + 32(p 2 /4 + l)/p 4 ] (= 0.013) 
and W = 2n^H = p, in good agreement with the heights and widths observed in Fig. |] and 
HI The zeros are expected to produce isolated peaks up to a value n = n c where their 
mean level spacing (measured in units of pit)) becomes comparable to the width W of 
the peaks. This indicates the onset of a resonance overlap regime. This condition gives 
n c fa log(i/27r) exp(27r) (~ 13100 for the value of t used in the numerical calculations). On 
the other hand, the contribution of the trivial zeros is small compared either to the complex 
zeros or to the pole due to their greater distance with respect to the imaginary axis. The 
curve with squares in Fig. [| is a plot of the r = 1 term of Eq. (|3.6| ) including all the singular 
points, compared to C{n) computed from Eq. flOP . The presence of negative resonances at 
the (rescaled) position of the lowest zeros was already noticed in the two-point correlation 
function of the Riemann zeros in Ref . fl20|) (see also the concluding remarks). 

Consider now the contributions of the singular points with r > 2, which are at the origin 
of the smaller oscillations in Fig. [5]. Their displacement in the complex plane as a function of 
r is given by Eq.( |3.7| ). In order to simplify the discussion, we consider explicitly the critical 
zeros n^ r , but the analysis can be extended to the other singular points. Two different 
effects on the location of these singular points with respect to the r = 1 distribution occur 
(see Eq. ( |3.7| )). On the one hand, as r increases their real part increases, and R,en^ r — > p as 
r — ► oo. This implies that the weight of their contribution to C(n) decreases with increasing 
r. The global factor (1 — r)/r 2 in Eq. fl3.6p goes in the same direction. On the other hand, 
their location is compressed towards the real axis as r increases, since Imn F = pt^/r. 
Finally, remember that the terms with r > 2 in ( |3.6| ) have the opposite sign with respect to 
the r = 1 terms. 

We can therefore summarize the contribution of the critical zeros in Eq. Q3.6| ) as follows. 
Each critical zero produces a series of peaks in the autocovariance. The main one has 
a negative sign and is located at n — pt^. Decorating the main resonance there is a 
set of smaller positive peaks located at n — pt^/2, pt^/3, . . . (with height diminishing as 
H — (r — l)(4vr 2 r 2 )- 1 log[l + 2((Re?v) 2 + l)/(Re?v) 4 ]). Therefore, at a given value of 
n, the main contributions to C(n) come, with decreasing importance, from zeros located 

at tn = n/p, 2n/p, 3n/p, The smaller oscillatory structure observed around a given n 

is thus due to zeros located r = 2,3,... times higher up the critical axis, thus producing 
an interference pattern to which a large number of nontrivial zeros are contributing with 
decreasing weight. The curve with diamonds in Fig. R is a plot of Eq.(B^) including up to 
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the r = 3 terms. We observe that including terms higher than r = 1 the agreement with the 
result obtained from Eq.( [3.1|) has greatly improved, and that the most important features 
are reproduced. Even peaks up to the third generation are easy to identify. For example, 
the small positive peak located at n « 18.4 corresponds to the third-order sub-resonance of 
the first Riemann zero pt\/3, with t\ ~ 14.13. The resonance formula Eq. ( |3.6[ ) thus allows 
for a simple and transparent interpretation of the structure of C(n). In contrast to the 
periodic orbit sum Eq.( |3.1| ), it allows to control the amount of complexity to be included 
hierarchically. 

Another interesting feature are the anomalously large peaks observed occasionally in 
C(n), as in part (b) of Fig. [|. As mentioned before, at the height t fixed by the actual 
numerical computations each complex zero of ((s) produces a peak in C(n) of height H « 
0.013 up to values n < n c . However, this is true on average, and statistical fluctuations of 
the spacings between consecutive zeros may produce interference effects between neighboring 
peaks. In some extreme cases, we may have two zeros which are very close to each other 
(on the scale of the width of the peaks). The occurrence of zeros whose distance is much 
smaller than the mean spacing is called Lehmer phenomenon, and is closely related to 
the Riemann hypothesis®. Whenever this unlikely phenomenon between two zeros occurs, 
their contributions to C(n) will add coherently and produce a peak twice larger than the 
contribution associated to an isolated zero. To illustrate this, we display in Fig. ^| the real 
function 

Z(t) =exp(i0(t)) C(l/2 + it) , (3.8) 

with 9{t) = Imlogr(l/4 + it/2) — (tlog7r)/2. The presence of two almost degenerate zeros 
at t — 7702/p = 1977.24 can be observed. This explains the large peak in C{n) of size 
H « 0.026 observed at n = 7702 in part (b) of Fig. |. 

The remarkable fact that, as for the optical image produced by two mirrors face to face, 
the low-lying zeros of the Riemann zeta function emerge as a decaying infinite sequence 
of replications when looking at the correlation of zeros located at an arbitrary height t on 
the critical axis appears more as a conspiracy than as a generic feature of a dynamical 
system, although this deserves further study. The origin of this singular property seems 
to be associated to three (apparently unrelated) peculiarities of the Riemann zeros which 
conspire to produce this effect: (i) the trace formula for their density is exact, (ii) the period 
of the periodic orbits, T p = logp, does not depend on energy (energy corresponds here to 
the location t on the critical axis), and (iii) the monodromy matrix has only one expanding 
eigenvalue. 

IV. SUMMARY AND CONCLUSIONS 

In the present paper autocovariances C(n) of two spacings between consecutive levels 
located n levels apart (n = 1, 2, . . .) are studied. Whereas the two-point correlation function 
i?2(e) can be viewed as the curvature of the number variance S 2 , C(n) represents (to within 
a sign) a discrete version of it. This fact implies that the effects present on the form factor 
(Fourier transform of R2) on the scale of the Heisenberg time T# are absent on C(n) (or 
on its Fourier transform). In particular, we find that the non-analytic oscillatory structure 
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of the RMT result of i?2(e), usually associated to off-diagonal contributions in semiclassical 
theories, are absent in C (n) which is, in the universal regime, a monotonic function. The off- 
diagonal semiclassical contributions are here expected to reproduce the smooth higher-order 
corrections written in Eq. fl2.3lD . 



We first derive expressions for the autocovariance within the RMT of Wigner-Dyson. 
They are approximate but accurate. In particular, it is shown that the autocovariances in 
RMT are negative, monotonic functions of n, and that for large n they vanish like — 1/ f3it 2 n 2 . 
The expressions found can be compared to those obtained using a similar but alternative 
route proposed in a previous paper. Both methods have their own advantages and draw- 
backs. The method adopted here allows to express systematically selected information on 
the form factor in terms of covariances. This is first illustrated for the case of quantum sys- 
tems which are classically diffusive. An explicit expression of the autocovariances is given 
which consists of a universal (random matrix) contribution plus a diffusive term which de- 
pends on the dimensionality d of the system. The latter, which appears to be identically 
zero for two-dimensional systems, dominates the large n behavior of the autocovariances 
and has opposite sign for d — 1 and d = 3. Large distant spacings are negatively correlated 
(anti-correlated) in one dimension and positively correlated in three dimensions. Finally, 
for any dimension if one neglects the finest time scale given by the elastic mean free time, 
autocovariances tend to zero in the large n limit. 

Section II. C (study of quantum systems which are classically chaotic) and section III 
(application to the zeros of the Riemann ^-function) constitute the main body of the paper. 
By using periodic orbit theory we arrive at an accurate formula for the autocovariances 
C(n) in terms of a diagonal approximation (Eq. ( [2.32] )). The important time scales are the 
Heisenberg time (which is taken as unity) and the period of a short periodic orbit r min . 
Based on this result, we then find an expression of the autocovariances in terms of a sum 
over derivatives of the usual (classical) dynamical zeta function Z(s) and another (classical) 
zeta function F(s) closely related to it (Eq.( |2.39D ). The interpretation of the singular points 
of Z(s) is well known (relaxation of an initial classical cloud, as given by the eigenvalues of 
the Perron- Frobenius operator). In contrast, we are not aware of a similar discussion of F(s). 
Its study remains an open problem, in particular its analytic structure and interpretation. 
Presumably that as far as autocovariances (and the two-point correlation function) are 
concerned, Z(s) gives the main behavior and the role of F(s) is limited to provide some 
additional finer structure. At any rate, it is the location in the complex plane of the singular 
points (zeros and poles) of Z(s) and F(s) that determine the structure of C{n) through a 
resonance- type formula (Eq. ( |2.43| )). The location of the classical singular points are rescaled 
by the Heisenberg time Tjj = hp. For ergodic systems, in the semiclassical limit Tjj — > oo, 
only the zero of Z(s) at the origin gives a non-negligible contribution and C(n) takes the 
universal form ( 2.14Q . At finite T# the non- universal singular points contribute significantly. 



The expression of C{n) in terms of classical zeta functions is easy to extend to the two 
point correlation function ^(e) = 1 — Y 2 (e). Using Eq. (|2.9| ) as starting point, and by similar 
steps as before we obtain, instead of Eq. fl2.31p , the result 

- 1 = £ E | det( 4-Di C08(2 " rT > ) + {R?MT{() + (4 ' 1) 

where we are assuming an average over a small energy window that kills some highly os- 
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cillatory frequencies in the sum over periodic orbits. i?^ MT is the usual two-point function 
of RMT. As mentioned in Section II. A, the random matrix correction to the diagonal sum 
involves a non-analytic oscillatory behavior contained in i?^ MT (e) (cf Eq. fl2.20D ), and which 
is absent in C(n). Due to this and contrary to the autocovariances of spacings, the diagonal 
approximation is a poorer approximation for R2. The term (/37r 2 e 2 )~ 1 plays here the role 
— C° m (n) is playing for the autocovariances, and compensates the smooth large-e decaying 
behavior of i?^ MT (e). The other important difference between ( |4.1[ ) and Eq. fl2.31| ) is the 
presence of the square of the sine function in the latter, which is due to the discrete nature 
of the curvature C(n). 

Exactly as for the autocovariances we may express the diagonal sum in Eq. ( [4.1| ) in terms 
of classical dynamical zeta functions 



*-» (£) _ 1 = _L Re |l l0! 



Z (ie/hp) 
F(ie/hp) 



(4.2) 



The relation connecting the diagonal part of R2 and Z(s) was found in Ref . (|T§D (where an 
extension to include the non-diagonal part was also given, see also Ref. (p4|) ) . Here we have 
generalized this by including a second dynamical zeta function, which we have explicitly 
computed in the case of the Riemann ^-function and turns out to bring non-negligible 
contributions. Its contribution was usually expressed as a sum over periodic orbits, and not 
written in terms of the more transparent approach used here. In contrast to Eq. ( |4.2| ), the 
existing relation between C(n) and the dynamical zeta functions Eq. (|2.39|) involves a sum 



over derivatives, whose first term coincides with Eq. fl4.2p . The summation leads to different 
shapes of the resonances of C(n) compared to those of Eq. fl4.2|) . 

These two quantities - the spacings autocovariance C(n) and the spectral two-point cor- 
relation i?2(e) - illustrate the manifestation of classical resonances in quantum correlations. 
This connexion seems however to be more general, as shown by a recent experiment© on 
microwave scattering in an open system. 

We finally study, in the framework developed in the previous section on chaotic systems, 
the autocovariances C (n) for the case of zeros on the critical line of the Riemann ^-function 
(section III). This is achieved by using the well known "dictionary" which translates quan- 
tum (eigenvalues of the Hamiltonian or evolution operator) and classical (periodic orbits, 
actions,...) information into the corresponding function quantities, namely imaginary part 
of zeros ("quantum") and prime numbers ("classical") (cf Table I). 

We first derive a periodic orbit sum (over primes and its repetitions) expression of the 
autocovariances C(n). After a short universal random matrix regime, C(n) is dominated 
with increasing n by the non-universal terms, encoded here in the prime numbers. We check 
the accuracy of the formula obtained by comparing to partial existing "data" of Odlyzko. 
The agreement is good. It is unclear whether the remaining differences are due to the 
approximations made in our derivation or to numerical inaccuracies. Probably to both. 
This deserves further study. 

Better insight can be gained by writing the autocovariances in terms of the (classical) 
dynamical zeta functions Z(s) and F(s). Extending the existing "dictionary" which connects 
dynamical properties and ^-function properties, one finds that the classical dynamical zeta 
functions Z(s) and F(s) are expressed in terms of the Riemann ^-function itself (Eqs.( fOl) 
and fl3.3|), and Table I). This remarkable connection between classical operators and the 
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function that determines the "quantum spectrum" leads to a resurgence phenomenon of 
the low lying Riemann zeros in the autocovariances computed at a height t » 1 on the 
critical line. This phenomenon is clearly displayed in a resonance formula obtained for C(n) 
(Eq. (|3.6| )) which is particularly appealing and which allows to interpret most of the features 
observed from data as predicted by the periodic orbit expression. 

The salient qualitative and quantitative features of the autocovariances C(n), mentioned 
in what follows, are easy to read off from the resonance formula. They have been checked 
by comparing to the results obtained from the periodic orbit summation formula, which 
is accurate but comparatively opaque. At a given height t on the critical line, which sets 
the scale p, C(n) starts as the random matrix prediction. Up to a critical value n c (oc 
logt) of n, the low lying zeros of ( produce (negative) well isolated constant-amplitude 
peaks centered at pt^, where t M is the imaginary part of the /i-th Riemann zero. One 
therefore neatly sees the low-lying zeros of ( appearing in the autocovariances computed at 
an arbitrary height on the critical line. The effect of the other dynamical zeta function, F(s), 
is to provide further structure, namely smaller (positive) peaks at values pt^/2, pt^/3, 
with decreasing importance. One therefore sees the low-lying zeros lurking as sub-resonant 
(positive) peaks. Beyond n c , the overlapping resonance regime sets in, and C(n) oscillates 
erratically. Occasionally, when two zeros of ( are very close (Lehmer phenomenon) at a 
height t, they give rise to a peak with double amplitude at n — pt. 

Let us finally conclude by noticing that in several respects it is easier to apply the 
formalism developed in Section II to the Riemann case, as we have done here, than to a real 
dynamical system. It may be worth to study in some detail some of the best known maps, 
like the cat or the baker map, from the present viewpoint. The physical significance of the 
dynamical zeta function F(s) introduced here may then become clearer. 
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TABLES 





Correspondences 
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Dp — > I lOg p 


pel 1UU 


lp ► iogp 


rescaled period 


— ► logp/ log(t/27r) 


Lyapounov exponent 




stability factor 


|det(M p r - 1)| 


dynamical zeta 


Z(s) ^ C -i(l_ s) 


functions 


^)-n~ 2 [c(r-r S )](- 1 )/'- 2 


"ergodic" zero of Z(s) 


70 -> pole of C(s) 


other zeros and poles 


— > zeros of £(s) 


of Z(s) and F{s) 





TABLE I. Correspondences for the "Riemann dynamics". The last three entries are discussed 
in the text. 
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FIG. 1. Spacing autocovariances for a three-dimensional diffusive system, with (3 = 2. The 
dashed curve is the random matrix term C rm ; the continuous curve includes the contribution C^a , 
Eq.(|2T25|). 
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FIG. 2. Spacing autocovariances C(n) in the range n < 20. Triangles: random matrix result, 
from Eq.( |2.16] ), Squares: data from Odlyzko computed from Riemann zeros in an interval starting 
around the 10 12 -th zero at t = 267653395648.8475 and containing 50000 zeros. Circles: theoretical 
curve, from Eq.( |3.1| ). Points are joined by solid lines to guide the eye. 



24 




25 



0.02 



C(n) 



0.00 



-0.02 



0.02 



0.00 



-0.02 
7600 

0.02 



0.00 



-0.02 
9800 




100 



200 



7700 



9900 



(a) 



300 



400 




7800 




n 



10000 



FIG. 4. Spacing autocovariances C(n) for Riemann zeros for three different ranges of n to 
illustrate several regimes and phenomena. Circles, from theory (Eq.( |3.1| ), see caption of Fig. |2|). 
The last part of Fig. 4(c) corresponds to Fig. |3[ 
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FIG. 5. Spacing autocovariances C(n) for the Riemann zeros. Enlarged view of the left part of 
Fig. 0(a). Circles: from periodic orbit theory, as in Figs ^, ||] and ||. Squares: from the resonance 
formula Eq.( |3,6| ) where only the r = 1 term has been kept. Diamonds: from Eq.( |3.6| ), including 
up to the r = 3 term. The two arrows in the lower part of the figure indicate the rescaled position 
pt^ of the first two complex zeros 1/2 + it a of C( s )- The eight in the upper part correspond to 
the first sub-resonances pt^/2 of the first eight zeros of £( s )- Some other (positive) peaks visible 
on the figure correspond to higher order sub-resonances. 
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FIG. 6. The function Z(t), Eq, (|3.8| ), showing two close-lying zeros of ((1/2 + it) around 
t = 1977.2. This almost degeneracy is at the origin of the large peak observed in C(n) at 
n = pt = 7702 in Fig.4(b). 
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